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Abstract. We give a direct proof of well-posedness of solutions to gen- 
eral selection-mutation and structured population models with measures 
as initial data. This is motivated by the fact that some stationary states 
of these models are measures and not functions, so the measures are 
a more natural space to study their dynamics. Our techniques are based 
on distances between measures appearing in optimal transport and com- 
mon arguments involving Picard iterations. These tools provide a sim- 
plification of previous approaches and are applicable or adaptable to a 
wide variety of models in population dynamics. 



1. Introduction 

Selection-mutation equations are models for structured populations with 
respect to continuous phenotypical evolutionary traits. They are usually 
written as equations for densities on the parameter space of phenotypes 
[3, 5, 19, 6], that is, they are usually formulated in spaces. However, for 
some of these models a more natural space to study the time evolution is 
the space of positive measures, since it has been proven in some cases [5, 13] 
that for a small mutation rate the steady states tend to concentrate in a 
Dirac mass at the evolutionarily stable strategy value. More generally, the 
need for a theory of well-posedness in measures for structured population 
models was mentioned in [15, 14, 28]. 

Some efforts in this direction have been directed at particular models 
in population dynamics: pure selection models for phenotypic traits in the 
space of measures were recently studied in [1, 11], while in [4] a particular 
case of a selection-mutation equation for a genetic trait in the space of 
measures is analyzed. More recently in [10] the author shows well-posedness 
and studies asymptotic behavior of a selection-mutation equation. 

Our aim here is to give a simple and general proof of well-posedness in 
the space of measures for a class of models that will include a wide range 
of selection-mutation models as well as many classical nonlinear structured 
population models [3, 14, 15, 22, 28, 27]. The basic model we consider is an 
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abstract Cauchy problem of the form 



dtu + V ■ {F{x)u) = N{t,u), 
u{0,x)=uo{x) (xGM'^). 



(la) 
(lb) 



where u = u{t, x) is the unknown, which depends on t > and x G M*^ 
(in any dimension d > 1). Any differential terms in the equation should be 
included in the term V • {F{x)u) (e.g., growth or aging terms); the term 
N(t, u) (a measure depending on x) may include nonlinear birth and death 
rates and selection- mutation interactions. Several examples in population 
dynamics are given in Section 3, where it is also shown how the abstract 
result may be applied to equations whose domain is not the whole space M'^. 

Our proof of well-posedness in the space of measures to (1) will be based 
on techniques stemming from mass transportation and semigroup theory. 
Similar ideas were already used for studying the mean-field limit of kinetic 
equations such as the classical Vlasov equation [16] and more recently in 
swarming models [8]. The main advantage of our approach is its simplicity, 
coming from the use of techniques already well developed in other fields, 
and its flexibility, which allows it to be adapted to a wide range of models 
in population dynamics. 

The well-posedness of measure solutions to some equations of the form 
(1) has recently been analysed by different although related techniques in 
[18, 17]. selection-mutation models were not included in their formula- 
tion but Sharpe-Lotka-McKendrick-type models (age-structured models) in 
which the boundary condition is introduced as a measure- valued right-hand 
side of the equation are treated, see subsection 3.3 for related results. 

A nonlinear semigroup approach using the splitting method for the trans- 
port V • {F{x)u) and the right-hand side N{t, u) terms was introduced in 
the more recent paper [9] to treat equations of the form (1). Here, we give a 
direct and simpler proof based on Picard iterations in the right metric space 
to conclude existence, uniqueness, and continuous dependence without re- 
sorting to the splitting method or nonlinear semigroup techniques. 

The organization of the paper is as follows: in Section 2 we show that the 
abstract Cauchy problem (1) is well posed for measures as initial data in the 
so-called bounded Lipschitz distance under reasonable Lipschitz conditions 
on F and N similar to the ones needed in [9]. Then, Section 3 is devoted to 
the application of these results to more explicit examples, mainly some non- 
linear selection-mutation models in subsections 3.1 and 3.2 where -F = 0, but 
also some mixed nonlinear structured population/selection- mutation mod- 
els in subsection 3.3 and pure structured population models in subsection 
3.4 where F ^ 0. These applications highlight the wide applicability of the 
abstract theorem in Section 2 for this type of models, setting a possible func- 
tional framework for stability and asymptotic convergence towards measure 
solutions [1, 5]. 



2.1. The bounded Lipschitz norm. Let us start with a quick summary 
of the definition and properties of the bounded Lipschitz norm, also called 
fiat metric [23, 25]. We denote by X(M'^) the set of Radon measures m M , 
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and consider the space of Lipschitz functions 1^^'°°(]R'^) endowed with the 
llV'lloo + Lip(?/'), with Lip{'tp) the Lipschitz constant of ip. 



Definition 2.1. Given a Radon measure fi G A4{W^) we define its bounded 
Lipschitz norm ||A*||_A4(]Rd), ||^|| when there is no ambiguity, by 



l^ll := sup 



dfi 



where C is the set of bounded functions : M'^ ^> M with Lip('0) < 1 and 
ll^lloo < 1, i.e., £ := {V € W^i'-(R'^) | HV'lli,^ < 1}. 

One sees from the definition that this is just the dual norm of W^'°°{W^), 
and that by duahty for any V G VF^'°°(]R'^), 



ipdfi< ||V'||i,oo||/^||. 

We remark that on probabihty measures this can also be defined as a dis- 
tance of Kantorovich-Rubinstein, or Wasserstein, type: when fi, v are prob- 
ability measures it holds that 

ll/i — z/|| = inf / min{|x — ?/| , 1} (i(T(x, y) , 

where S is the set of transference plans between fi and z^, that is, probability 
measures on M*^ x with marginals fi and u [25]. 

We will need to use the following simple result on this distance: 

Lemma 2.2. If b e and fi G M{M.'^), then bfi G A^(IR'^) and 

\\bfJ-\\ < ||&||l,oo||/W||. 



Proof. It is clear that bfi G A4{M. ), since 6 is a bounded continuous function. 
We integrate against tp G W^'°° with HV'lli oo < 1 to find 



bip dfi 



< II^V^II^ooll/zll (2) 



On the other hand, we have 



so 



oo < ll&lloo II V^lloo, 



Lip(6V) < Lip(6)||V'||oo + ||6|ULip(V'), 



|i,oo = ||6V'l|oo + Lip(6V^) < ||6||oo(||V'l|oo + Lip(V')) +Lip(6)||V^|U 
< ||6||oo + Lip(6) = ||6||i,oo. 
Plugging this into eq. (2) finishes the proof. □ 

In the rest of this paper, we will work with measure solutions to some 
evolution partial differential equations and therefore, we will work with the 
space of bounded continuous curves on the set of measures i?C(/; 7W(M'^)) 
depending on t G / denoting the time variable, with / = [0, T] for some 
T > or / = [0, oo). The continuity of the curves of measures t i— ?> fi{t) is 
always understood to be with respect to the bounded Lipschitz norm. We 
warn the reader that elements in BC{I; A4(M'^)) will often be denoted as if 
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they were absolutely continuous densities with respect to Lebesgue with the 
form dfi{t) = u(t, x) dx for the sake of simplicity. 

The standard total variation norm for measures will be denoted by || • ||tv- 
We remark the natural necessity of the bounded Lipschitz distance (or simi- 
lar distances between measures) to work with transport evolution equations, 
as opposed to the total variation norm. In fact, take any injective contin- 
uous path X : [0, T] — > M*^ and take the curve of measures fi defined by 
t — > 6x(^t)- It is easy to check that ^ belongs to SC([0, T]; X(]R'^)) while 
\\n{t) - ^^{s)\\TY = 2 for all < t < s < T. 

Although all models in population dynamics study the evolution of posi- 
tive measures (number density of individuals with respect to some variables) , 
let us mention that we need to use the bounded Lipschitz norm and not other 
optimal transport distances since the total mass (total variation) of measure 
solutions will typically not be preserved in time. We will denote by Bbl{R), 
resp. BTYiR), the ball of radius R centered at in the Bounded Lipschitz 
norm and in the total variation norm resp. Finally, let us point out that 
balls Bt^y{R) in Ai{M.'^) with respect to the total variation norm are closed 
in the bounded Lipschitz norm by simple weak convergence arguments. 

2.2. An abstract result. We consider the abstract evolution equation for 
measures given in (1), which we recall here: 



Here u = u{t, x) is the unknown, which depends on t > and x G M.'^, under 
the following hypotheses on F, N and uq: 

(HI) uo e M{R'^). 

(H2) F : M'' — )• M'^ is a bounded Lipschitz map. 

(H3) N : [0, +oo) x MiW'^) A4(M"') is a continuous function both in 
t and u (considering here the topologies induced by the bounded 
Lipschitz norm on A4{M.'^) and the usual topology on [0, +oo).) 

(H4) N is locally Lipschitz in its second variable, i.e., for every bounded 
set K C [0, +oo) X M{W^) there exists = Ln{K) > such that 

\\N{t,ui) - N{t,U2)\\ < Ln\\ui-U2\\ V {t,ui),{t,U2) G K. 

(H5) N carries bounded sets in the total variation norm to bounded sets 
in the total variation norm: for each R > there exists Cr > 
such that ||A^(t,ii)||TV < Cr for all t > and u G M{W^) with 
||u||tv ^ R- 

Definition 2.3. Assume Hypotheses (H1)-(H5), and take T G (0, +oo]. We 
say u G BC{[0,T); Ai{W^)) is a solution of equation (1) on [0,T) with initial 
condition uq when, for every (p G Cq°([0, T) x R'^), 



dtu + V^. • {F{x)u) = N{t, u) 
u{0,x)=uo{x) (xGR^^). 




i 



■T 




F{x)V(j){t, x) u{t, x) dx dt 
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Theorem 2.4 (Well-posedness of the abstract equation). Assume Hypothe- 
ses (H1)-(H5). There exists a maximal time T > such that there is a 
unique solution u € B{^,T); M.{W^)) of equation (1). In addition: 

(i) Either T = +oo or \\m.t^T ||'u(*)||tv = +oo. 

(ii) This solution depends continuously on the initial condition in the 
hounded Lipschitz norm: take two solutions ui,U2 of equation (1) 
on [0,T) with initial conditions u^, respectively. Assume also 
that 

\\uiit)l\\u2m<R (tG[0,r)), 

and take to he the Lipschitz constant of N with respect to the 
second variable on the set [0, T) x Bbl{R) ^ [0, +cxd) x M{W^). 
Then, 

||ni(t)-n2(t)|| <e(^^+^^)*||n?-nO|| (tG [0,r)). 

Remark 2.5. If there is no drift term present in equation (1) (this is, F = 0) 
and C M'^ is an open set, then the resuh holds if one changes M{W^) 
by The modifications needed in the proof below are straightforward 

and we omit them. 

Proof of Theorem 2.4. Define Xt : ^ M!^ as the fiow at time t of the 
characteristic equations 

By standard arguments in the theory of ordinary differential equations we 
have that Lxt, the Lipschitz constant of the fiow Xt at time t, satisfies 

Lx, < e^^' (t > 0), 

where Lp is the Lipschitz constant of -F. 

We will prove existence by a fixed point argument in the set 

Mt ■■= {u G BC{[0,T];Bry{R)) \ u{0) = u°}, 

of radius R := 2||no||TV- We choose 

T := min{||uO||Tv/CH,, l/Lp, l/(3Ljv)}, (3) 

where Cr and Ljy are given in Hypotheses (H1)-(H5) with K = [0,T) x 
Bbl{R)- We endow Mt with the standard norm 

|||n||| := sup \\u{t) II . 

tG[0,T] 

As noticed before ^AT is a closed subset of the space i?C([0, T]; i?Tv(-^?))- 
Hence Mt is complete and we may apply the Banach fixed point theorem 
in it. Define the map F : ^AT Mt as 

V{u){t) := Xt^u"" + [' Xt^s#N{s,u{s))ds, (4) 
Jo 

with ^ denoting the push- forward of a measure through a map. It is easy 
to check that a fixed point of this map is in fact a solution to equation (1). 
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Step 1: r is well-defined.- Let us first show that t i— )■ Xt^u^ and 
s I-)- Xt-s#N{s,u{s)) (for t € [0,T] fixed) are continuous maps. For the 
first one, taking any test function (/> G £ and any t,T £ [0,T], 

/ • - dx= I (HMx)) - Cl)i^rix)))u\x) dx 

< [ \Xt{x) - Xrix)\ \u^\{x)dx 

< |t - t| ||-F||oo||^i°||TV, 

which shows continuity. As for the second one, we take (j) as before, fix 
t € [0,T] and take any r, s € [0,t]. Denoting N{s,u{s)) as Ns for short, we 
have 

/ • {Xt^s#Ns - Xt-r#Nr) dx 

= I (t>- {Xt.s#Ns - Xt^s^Nr) dx+ [ Ct>- {Xt-s#Nr " Xt-r#Nr) 



dx 



:4>oXt^s){Ns-Nr)dx+ / {(f>{Xt-s{x)) - (l){Xt-rix))) Nr{x) dx 

< Lx,_J\Ns - NrW + \t- s\ ||F||oo||iVr(x)||TV 

< e(*-^)^^ \\Ns -Nr\\ + Cr||F||oo|t - s\ . 

This proves continuity, as s i— )> N{s,u{s)) is continuous due to (H3). Hence, 
the integral in (4) makes sense, r(n) is continuous from [0, T] to A4{M.'^) in 
the bounded Lipschitz norm, and we only need to see that its image is inside 
Btv{R): 

||r(n)(t)||TV < \\Xt#u%^+ f \\Xt-s#Niu{s))\\^y ds 

Jo 

< ||nO|| + / ||iVKs))||Tv ds 







I TV 

< \\u^\L,, + CrT < 2|U°|L,, = r. 



ITV ' ^« - II IITV 

Step 2: r is contractive.- Take u,v Mt- Using similar arguments 
we estimate 

\\r{u){t)-T(v){t)\\ < r \\Xt^s#N{s,u{s)) - Xt-s#N{s,v{s))\\ds 
Jo 

< I Lx,_AN{s,u{s))-N{s,v{s))\\ds 
Jo 

< c^^'^Ln / Ms) -v{s)\\ds. 

Jo 

By taking the maximum over t G [0, T] this implies 

|||r(u) - r(7;)||| < c^p^LnT \\u - v\\ < L ju - v\l 

for some L < 1, due to the choice of T made in (3). An application of the 
Banach fixed point theorem together with usual arguments on the extension 
of solutions finishes the proof of point i) of the theorem. 



MEASURE SOLUTIONS FOR SOME MODELS IN POPULATION DYNAMICS 7 



Step 3: Continuous dependence.- We estimate the difference of the 
two solutions as follows: 

\\u{t) - vit)\\ 

<\\Xt#u^ -Xt#v'^\\+ [ \\Xt^s#N{s,u{s)) - Xt-s#N{s,v{smds 

Jo 

<LxJk°-^^°||+ / Lx,_M{s,u{s)) - Nis,v{smds 
Jo 

<e^^^\\u° -v°\\+Ln [ e^^^^-'^\\u{s)-v{s)\\ds. 
Jo 

Gronwall's Lemma then implies the result. □ 

Remark 2.6. Theorem 2.4 is a generalization of ideas in the theory of linear 
evolution semigroups [21, 2], since equation (1) is the sum of a linear term, 
and a locally Lipschitz perturbation. However, a small modification of the 
argument is needed: this comes from the fact that one cannot work in the 
dual space [M^^'°°(M'^)]*. Actually, the proof above shows that by restricting 
to measures in i?Tv(-R) we are able to prove the continuity of the transport 
semigroup. This continuity is not evident in Bbl{R)- 

Finally, we point out that we are usually interested in positive measures 
as initial condition, even if Theorem 2.4 does not require positivity. It is 
most often the case in models that the quantity under study is the density of 
a given population, which is intrinsically positive. Hence, for most models 
of interest, positivity is preserved in time (see for example Lemma 3.2). 

3. Application to particular models 

In this section we will apply Theorem 2.4 to show well-posedness of four 
particular models in population dynamics. The first one is a simple selection- 
mutation equation for a phenotypic variable inspired by the "continuum of 
alleles model" introduced by Crow and Kimura (see [12] and also [4]) in the 
field of population genetics in order to explain the maintenance of genetic 
variation due to the balance effect of selection and mutation. The second 
one, introduced in [6], is a modification of the first one in which it is assumed 
that the nonlinear term modelling the competition between individuals for 
resources is infinite-dimensional. The third model we consider was studied 
in [7] and it is a selection-mutation model for an age-structured population. 
The last model we present is an age- and size-structured model that was 
introduced in [26]. 

3.1. A simple selection- mutation equation. Let us consider the follow- 
ing selection-mutation equation: 

— (t, x) = {1 — e)b{x)u{t, x) — m{x, P{t))u{t, x) 

+ e I b{yh{x, y)u{t, y) dy := N{u{t, ■)){x) (5) 
Jq 

for the density u{t, x) of individuals at time t > with respect to an (ab- 
stract) evolutionary variable x in an open set Q, C R'^. P(t) denotes the 
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total population at time t 

P{t) := / u{t,x)dx 

and m is the trait-specific death rate which depends in an increasing way 
on the total population P{t) at time t. The inflow of non-mutant newborns 
will be given by (1 — e)6(x)n(t, x) where b{x) is the trait-specific fertihty and 
e stands for the probability of mutation. The inflow of mutant newborns 
will be given by the integral operator e b{y)'y{x, y)u{t, y) dy where 7(x, y) 
is the density of probability that the trait of the mutant offspring of an 
individual with trait y is x. 

We may apply Theorem 2.4 to equation (5) under the following conditions: 

Theorem 3.1. Assume (HI) and also that b, m and 7 satisfy the following: 

(i) 6 : $7 ^ M is in W'^'°° (i.e., it is bounded and Lipschitz). 

(ii) m : X M M satisfies that for each p G M, m{-,p) G W^'°° , and 
for each R > there exists Lm > such that 

\\m{-,pi)-m{-,p2)\\i^^<Lm\pi-p2\ for allpi,p2 €[-R,R]. (6) 

(iii) For each y G ^l, 7(-,y) is a positive probability measure on Q, and 
there exists > such that 

Il7(-,?/) - 7(-,^)ll < ^7ly - ^1 forally,z£n. (7) 

Then the operator N in equation (5) satisfies the hypotheses of Theorem 2.4. 
Consequently, equation (5) is well-posed in the sense of Theorem 2.4. 

Proof. We need to check that assumptions (H4) and (H5) are satisfied, as 
(HI) is included in the statement, (H2) is trivial here since F = 0, and 
(H3) is a consequence of (H4) since the operator in equation (5) does not 
depend on time. We point out that due to Remark 2.5 we may use Theorem 
2.4 in 7W(0), as we have no drift term here (-F = 0). 

For (H4) we need to prove that given R > there exists a constant L > 
such that 

\\N{n)-N{iy)\\ < 
for all /U, G with ||^||, < R. For the first term in (5), 

Wbfi - bi^W = \\b{iy - /i)|| < ||6||i^oo||i^ - (8) 

For the second. 



< 



\{m{-,P{fi))-m{-,P{u)))fi\\ + \\m{-,P{u)){fi - u) 



< ||m(-,P(/u)) - m(-,P(i/))||i,oo||Ai|| + ||m(-,P(z^))||i,oo||^ - v\\ 

<K^P{|I)-P{v)\M + C\\^l-v\\ < \\fi-iy\\iC + LM\), (9) 

where C is a constant such that ||?7i(-,p)||i^oo ^ C for all p G [—R,R\ finite 
due to (6). Finally, in order to estimate the third term we notice that, for 
all £ C, 

-i{x,-)'4){x)dx <C, (10) 
1,00 
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for some C > 0. Indeed, / 7(x, y)^{x) dx is uniformly bomided for y G 
due to the fact that 7(-, y) is a probabihty measure, and it is also Lipschitz 
in y since 



ilix,y) - -i{x,z))'4}{x)dx 



< ||7(-,y)-7(-,^)ll <^7ly-^l 



for all y,z il,, due to (7). Hence, (10) holds and we can estimate the third 
term in (5) by integrating against a function ■0 S £: 



bivhix, y)iKy) - ^{y))i^{x) dy dx 



iKv) - '^iy))Ky) / l{x,y)ilj{x)dxdy 



< 11^* - '^ll ll^lll,. 



7(2;, •)0(x) dx 



<C\\l^-v\\ \\b\\. 



1,00 



Putting together (8), (9), and (11) we conclude that (H4) holds. Finally, 
(H5) is easily seen to hold using that b is bounded and 7(-, y) is a probability 
measure. □ 

In the general abstract theorem we do not show conservation of positiv- 
ity for solutions. Since L^{Q) is dense in in the bounded Lipschitz 

distance, this is a straightforward consequence of the result of conservation 
of positivity of solutions, which is already available for all of the models 
mentioned here. We show positivity of solutions for this model for the sake 
of completeness, but this will be skipped for the rest of the models of the 
paper, to which analogous arguments are applicable. 

Lemma 3.2. Under the hypotheses of Theorem 3.1, every solution of (5) 
with positive initial condition uq is positive. 

Proof. We begin by showing positivity of local solutions of (5) in L}-{W^). 
The initial value problem can be written as 



du 
dt 



Au + f{u) 

m(0) = Uq 



(12) 



where 



Au{t,x) := {I - e)b{x)u{t,x) + e l b{y)-f{x,y)u{t,y) dy, 

Jn 

and f{u){t,x) := —m{x,P{t))u(t,x). The operator A is the generator of a 
positive semigroup T(t). Let A be a constant bigger than the bound of m. 
If we add and subtract Am to (12) we get 



du 
'dt 



{A-Xl)u + f{u{t)) + Xu{t), 



The mild solutions of this new initial value problem, and therefore also those 
of problem (12), can be constructed by iterating a suitable variation of con- 
stants formula [21, 2]. More precisely, they are limits of the sequence {zn)n>o 
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of functions defined on [0, tmax) for some tmax > 0, recursively defined by 
the formula 

{t)=T{t)zQ+ / T{t-s)[f{Zn{s))+\Zn{s))ds, 

Jo 

where T{t) is the semigroup generated by the operator A — XI, that is 
f{t) = e-^^T{t). 

Since zq is positive, the semigroup T{t) is positive and since A is larger 
than the bound of m, we obtain that zi is positive. By induction over n 
we have that {zn)n>o is positive. Finally, since the cone of the positive 
functions of is closed, we obtain that z{t) is positive. Positivity of local 
solutions implies positivity of global solutions by a standard connectedness 
argument. Finally, using Theorem 2.4, the density of L}{Q) in in the 

bounded Lipschitz distance gives us conservation of positivity in the space 
of measures. □ 

3.2. A selection-mutation model with infinite-dimensional environ- 
ment. Another example of a selection- mutation equation is 

^{t,x) = ({l- e)h{x) - do{x) - [ d{x,y)u{t,y) dy)u(t, x) 

^ ^ (13) 

+ e b{y)-f{x,y)u{t,y)dy =: N{u{t,-)){x). 



Jn 

for the density of individuals u{t, x) with respect to an (abstract) evo- 
lutionary trait x € 0. The difference with (5) is that here the trait- 
specific per capita death rate is given by the sum of the terms do{x) and 

d{x, y)u{t, y) dy. The latter one models the interaction between individu- 
als through competition for resources, and is the only nonlinear term in the 
equation (whose nonlinearity in this case is infinite dimensional). 

This model was presented in [6], where the authors prove existence of 
steady states and also that their asymptotic profile when the mutation rate 
e ^ is a Cauchy distribution. Our well-posedness result in the space of 
measures for equation (13) is the following: 

Theorem 3.3. Assume (HI), points (i) and (in) in Theorem 3.1, and also 
that do and d are nonnegative functions satisfying 

do:n^RisinW^'°°{n), (14) 

andd-Mxn^R with d € W'^''^{n; W'^''^{n)); that is, there exists L>0 
such that 

N(^r)lliyi,oo(Q) < L for all xe^. (15) 
\\d{x, ■) — d{z, ■)\\yyi,oo^Q'j < L\x — z\ for all x , z £ ft . (16) 

Then the operator N in equation (13) satisfies the hypotheses of Theorem 
2.4. Consequently, equation (13) is well-posed in the sense of Theorem 2.4. 

Proof. As remarked in the proof of Theorem 3.1, we only need to check (H4) 
and (H5). As the other terms have the same form as the terms in (5) since 
(14) is satisfied, we only need to check (H4) and (H5) for the term which 
involves d. 
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First we notice that due to (15)-(16) the term d{x , y)u{y) dy is in 
for any u G M{^1), as for all x G 0, 



d{x,y)u{y)dy 
and for any x,z G 0, 



< 



\u\\ < L\\u\\. 



{d{x,y) - d{z,y))u{y)dy 



< 



Actually, we have proved that for any w G 



d{-,y)w{y)dy 



< L\\w\ 



lull < L\x — z\ \\u\ 



(17) 



In order to prove (H4) for the term involving d, take two measures u, v in 
M{n). Then, 



u d{-,y)u{y)dy -V I d{-,y)v{y)dy 



< 



(u-v) / d{-,y)u{y)dy 



+ 



< \\u — v\ 



n 



d{-,y)u{y)dy 



n 
+ \\v 



V I d{-,y){u{y) -v{y))dy 

d{-:V){u{y) -v{y))dy 



1,CXD 



l,oo 



< Llln — i)||||n|| + L||f||||n — 



where we used (17) for the last step. This proves (H4). On the other hand, 
(H5) is easily proved since, in particular, \j^d{x,y)u{y)dy\ < L||n||TV for 
ah rc G r2. □ 

3.3. A selection-mutation model with age structure. Let us consider 
the following equation 



du du 

— {t,a,x) + —{t,a,x) 

ot da 

u{t, 0, x) 



— m{a, X, P{t), Q{t))u{t, a, x) 
il-e 



+ e 



oo roc 



n(0, a, x) = no(a, x) 



(18a) 

b{a,x)u{t,a,x) da (18b) 

7{x,y)b{a,y)u{t,a,y)dady (18c) 

(18d) 



where u{t, a, x) is the density of individuals with age o > and matu- 
ration age X > (the evolutionary variable) at time t. P{t) and Q{t) 
denote, respectively, the total population of juveniles and adults, that is 
P{t) = u{t,a,x) dadx, Q{t) = u{t,a,x) dadx, ni is the mor- 

tality rate, b is the fertility rate and 'j{x, y) is the probability density that 
the maturation age of the mutant offspring of an individual with matura- 
tion age y is x. As in the previous examples, e stands for the probability of 
mutation. 

This model is a slightly modified version of the one studied in [7], where 
the only difference is in the term of inflow of newborns. The difference 
of this model with the ones in the previous sections is that here, fixing 
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the evolutionary variable, we still have an infinite-dimensional model, more 
precisely, an age-structured population model. In [7] well-posedness of the 
model was proved in the Banach space L^{M.^) (denoting = [0,+oo) x 
[0,+oo)), and also the existence of steady states. In order to show well- 
posedness in the space of measures we rewrite (18) as follows: 

dtu + dau = Ni{u) (19a) 
u(t,0,x) = n2iu) (19b) 
u{0,a,x) = uo{a,x) (19c) 

where we call Q := and define, for u E A4{Q), 
Ni{u) := — m{a,x, P{u),Q{u))u, 



n2{u){x) := {1 — e) j b{a,x)u{a,x) da 

) 

7(x, y)b{a, y)u{t, a, y) da dy. 



Jo. Jo Jo 



The model may be rewritten in the form (1) by extending it to an equation 
on M^, with an additional independent term. Let us be precise about the 
intended solutions: 

Definition 3.4. Take T G [0,+oo]. We say a continuous function u : 
[0,T) — )• A4{n) is a solution of equation (19) on [0,T) with initial condition 
uq when, for every (j) G C^{[0,T) x Q), 

udt(j) dx da dt — / (j){0,x,a)uQ{x,a) dx da 
Jn 

T r j-T i-oo 

uda4>dxdadt— j j n2{u{t))(j){t,x^d) dx dt 

T 

Ni{u{t))(j)dxdadt. 

(When the variables oi (j) or u are not specified, it is understood that they 
are {t, a, x)). 

We now take a suitable extension of the functions m, b and 7 to all of 
(for definiteness, by mirror symmetry first in x and then in a) and consider 
the following equation, posed in the whole set of (a, x) G M^: 

du du 

^ + ^ = A^i(n) + n2{u)6a=o, (20a) 

u(0) = Uq. (20b) 

Equation (20) is of the form (1). Now, observe that a solution of (20), in the 
sense of Definition 2.3, is also a solution to (19) in the sense of Definition 
3.4 when restricted to M^, provided it is zero on the set Hence, we 

just need to give conditions on m, b and 7 so that (20) satisfies Hypotheses 
(H1)-(H5) and its solutions are supported on M^. 

Theorem 3.5. We assume the following: 
(i) b G W^'°°{Q,), and it is nonnegative. 
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(ii) m: QxM.xM^M.isa nonnegative function satisfying a condition 
similar to the one in Theorem 3.1; for each g € M, m{-,p,q) G 
W^'°°, and for each R > there exists L„i > such that 

\\'m{-,pi,qi) - m(-,p2,g2)|li,oo ^ Lm{\pi -P2\ + \qi - 92I) 

for all pi,p2,qi,q2 G [-R,R]. 

(iii) For each y G R, 7(-,y) is a positive probability measure on M and 
there exists > such that 

\h{-,y) -7(-,^)ll < ^7!^-^! /o'^ aUy,z e M. 

Then the initial value problems (20) and (18) are well-posed in the sense of 
Theorem 2.4. 

Proof. As remarked in the proof of Theorem 3.1, we only need to prove (H4) 
and (H5), and we can do it separately for each term. 

The term Ni{u) can be treated in a similar way to the term m{x, P)u in 
(5), and we omit the details. For the term n2{u)6a=o we have 

||(n2(M) - n2{v))6a=o\\MiR2) = \\n2iu) - n2(w)||_^(iK) , 

where the norm on A^(]R) and A4{M?) is the bounded Lipschitz norm. The 
term in n2{u) which involves 7 is of a similar form to the one in (5) and 
can be treated analogously. For the other term, taking any test function 
cj) G W^'°°{R) with ||0||i,oo < 1, 

(j){x)b{a,x){u{a,x) — v{a,x)) dadx 

< \\u - v\\ ||</'6||l,oo < \\u - v\\ ||6||l,oo, 

which shows that 



b{a, •){u{a, •) — v{a, •)) da 



< \\u — v\ 



1 1,00 J 



hence proving (H4) for this term. Condition (H5) for this term is easily seen 
to hold by using that b is bounded. 

The above allows us to apply Theorem 2.4 to equation (20). We deduce 
that the problem (20) is well-posed, and we only have to show that its 
solutions have support on M^, so that they are also solutions to (18). In 
order to do this, we take any time T > and any test function (pT G C°°(M^) 
with compact support on \ and consider cp : [0, T] x M x M — )• M to be 
the solution to 

^ - ^ = on (0, T) X M X M, (21) 
dt da ^ ' ^ ^ 

with (j){T,a,x) = (pTia^x) for a,x G M. Then, 

— / f (j)udadx= f I 4){—daU + Ni{u) + n2{u)5a=Q) dadx 

+ / / uda(t>dadx= / / (l){Ni{u) + n2{u)6a=o) dadx = 0, 
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since for all t € [0, T], (j){t, •, •) has support contained in M? \ M^, which can 
easily be seen since the solution to (21) is explicit. Hence, since the initial 
data is supported in M^, then 

(l)T{a,x)u{T,a,x) dadx = / / (j){0,a,x)uo{a,x) dadx = , 



and since (px was arbitrary we deduce that at time T, u has support con- 
tained in M2_. □ 

3.4. An age-size structured model. Let us consider the following age- 
size structured model from [26]: 

du du d 

+ ^ + ■Q^i9{x)u) = -m(a, x, P{t))u (22a) 
for a G (0,ai),rE G (xo,xi), t > (22b) 

u{t,0,x) = / / f3{a,x,x)u{t,a,x) dx da (22c) 
Jo Jxo 

for X G (xo,xi), t > (22d) 

u{0,a,x) = uo{a,x) for a G (0, ai), x G {xq,xi), (22e) 

where u = u{t, a, x) denotes the density of individuals with age a, with 
< o < ai < oo and size x with < xq < x < xi < oo. Size increases 
with time, in the same way for all individuals of the population, and the 
growth rate is given by the function g{x) which is assumed not to depend 
on environmental factors. Moreover it satisfies ^(x) > and g(xo) = 0. 
m denotes the mortality rate and /3(a, x,x) denotes the average number 
of offspring of size x produced per unit of time by an individual of age a 
and size x and P{t) = J^^ J^^ u{t,a,x) dx da. Here, we denote by = 
(0, ai) X (xq, xi) the domain of definition of the equation. 

Many versions of the model (22), both linear and nonlinear, have been 
studied, for instance in [20] and also in [24] where a more general nonlinear 
model containing an arbitrary number of structured variables is considered. 
The usual space to study these models is L^(0). By using essentially the 
same ingredients as in the previous subsection, one can prove the following 
theorem that we state without proof. 

Theorem 3.6. We assume the following: 

(i) m : $7 X R — )• M is a nonnegative function satisfying a condition 
similar to the one in Theorem 3.1: for each p G M, m{-,p) G W'^''^ , 
and for each R > there exists Lm > such that 

\\m{-,pi) - m(-,p2)|li,oo ^ Lm\pi -P2\ 

for allp\,p2 G [—R,R\. 

(ii) g G iy^'°°([xo,xi]) with ^(0) = and g{xi) > 0. 

(iii) The map /3 : — )• Ai{[xo, xi]) assigns (a, x) i-)- /3(a, x, •) and verifies 
thatW^'°°{il:,A4{[xo,xi])), i.e, it is hounded and there exists > 
such that 



(ai,xi,-) - /3(a2,£2, •)!! < Lis{\ai -02! + |xi -X2I 
for all (ai, xi), (02, X2) £ ^- 
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Then the initial boundary value problem to (22) is well-posed in the sense of 
Remark 3.7. 

Remark 3.7. Let us mention that the extension outside the reahstic domain 
Q to of the model ingredients m, /3, and g while meeting the conditions in 
Theorem 2.4 may be done in many different ways. Once one has an extended 
equation in M^, Theorem 2.4 applies, and all solutions to the extended equa- 
tions lead to the same solution once restricted to Q. This is due to the fact 
that the characteristics associated to the transport field {l,g{x)) for the age 
and size variables (a, x) are not incoming at the boundaries: a = ai, x = xq, 
and X = xi. A similar argument as in the proof of Theorem 3.5 shows that 
if the solutions for these extended systems are zero initially in the set of 
a < 0, then they remain so for all times. 
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